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Ecuaciones Diferenciales Ordinarias de 
Orden Superior 


1. Conceptos básicos 


En este tema estudiaremos ecuaciones diferenciales lineales de orden mayor que uno. En 
concreto, si / C ¡R es un intervalo de la recta real y a;, b: 1] —> R son funciones reales, 
para 0 <i< mm, con a, (1) 4 0 para todo x € [, nos centraremos en aquellas ecuaciones 
diferenciales que se puedan escribir de la siguiente forma: 


an (2) ul (2) +---+a2 (2) u” (2) +01 (2) u' (1) + 09 (2) u (1) =b(x) Vel. (P) 


A las funciones a; se les denomina coeficientes de la ecuación diferencial lineal. En 
el caso en el que a; (1) = a; € R para todo 0 < ¿ < n se dice que la ecuación diferencial es 
de coeficientes constantes. Cuando la función b = 0 se dice que la ecuación diferencial es 
homogénea, es decir, 


an (au (+ +0 (ou (0) +a (2) u (2) +0 (x)u(1)=0 Vel. (Pp) 


En lo que sigue, llamaremos (P) a la ecuación diferencial lineal no homogénea y (Py) a la 
ecuación diferencial homogénea asociada. 

En primer lugar analizaremos si una ecuación diferencial de este tipo tiene solución 
así como cuantas soluciones puede presentar. Para ello introducimos algunos conceptos y 
resultados que nos ayudarán a entender la naturaleza de las soluciones de (P). 


Teorema 1.1 Si ag, a1,...,a,,b son funciones continuas en I, con a, (1) 4 0 para todo 
rel ysizo € Í, entonces existe una única solución, u, del problema de valor inicial, 
PNL: 

Í An (2)UM (1) +---+a1 (2) u' (2) +00 (2) u(x) =b(x) 


u (2o) = 0%, ul (xp) =07,... UD (xp) = An-1 
con QA, Q1,...,An-1 ER. 
Definición 1.2 Sea TI C R un intervalo y sean fi, f2,..., fm : Il > R m funciones. Se 
dice que fi, f2,..., fm son linealmente dependientes en Í si 3C,,C2,...,Cm ER, no todos 


nulos, tal que 
cf (2) + cf 1(1)+-:+Cmfmíz) =0 Vx € Il. 


Se dice que fi, f2,..., fm son linealmente independientes cuando no son linealmente 
dependientes. 
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Puesto que las ecuaciones diferenciales (P) y (Py) son lineales, se deducen las siguientes 
propiedades: 


Proposición 1.3 La combinación lineal de dos soluciones de (Py) también es solución de 


(Pr). 


Demostración. 
Sean uz y us dos soluciones de (Py), sean a, PB ER y sea w = qu; + Bus. Entonces, 


= a | an (1) a (1)+---+a1 (2) u, (2) + a0 (2) uy ()) 
+8 (a, (x) y (1) +---+a1 (2) ua (2) + ap (2) ua (1)) =0. 
n 


Proposición 1.4 Si u, y uz son dos soluciones de (P), entonces w = us — uz es solución 
de (REL 


Demostración. 


(u, (2) — us (2) + 00 (2) (us (2) — ue (9) 
— Ue (2)) + do (2) (us (1) — us (2) 


Proposición 1.5 Siup es una solución de (P), entonces cualquier solución u de (P) verifica 
que u= UP + Uy, donde uy es una solución de (Pp). 


Demostración. Sea u una solución de (P) y sea uy =U— Up, entonces, cOmO U y UP SON 
soluciones de (P), por la Proposición 1.4 se tiene que uy es solución de (Py), y por tanto u 
se puede escribir como suma de up más una solución de (P,/). Ml 


Observación 1.6 Esta proposición proporciona un método para calcular todas las soluciones 
de (P) a partir de una solución particular (esto es, una función concreta que verifica (P)) 
y de todas las soluciones de (Py). Por tanto, en lo que sigue estudiaremos la estructura de 
las soluciones de (Py) y posteriormente describiremos métodos para encontrar una solución 
particular de (P). 
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Definición 1.7 Seal] CR un intervalo y sean fi, f2,..., fm: TI >R m funciones que poseen, 
al menos m—1 derivadas en I. Se define el Wronskiano de este conjunto de funciones como 
la función: 


O A 
tas PO A 
[AD () Pa) 2 Da) 
Proposición 1.8 Sea Il CR un intervalo y sean fi, f2,..., fm : TI >R om funciones que 


poseen al menos m — 1 derivadas en I. Entonces, 


fifa, Ffmy es linealmente independiente => 110€ 1tq.W (fi, f2,..., fm) (to) 40. 


Demostración. 
>>) Supongamos que [f,, fa,..., fm) es linealmente independiente. 
Razonando por el método de reducción al absurdo, supongamos que W (fi, fa,..., fm) (1) = 


0 Vx € 1. Puesto que el determinante que define el Wronskiano es cero para todo x € l, 
la última columna se puede escribir como combinación lineal del resto de columnas, lo que 
implica, en particular, que 3C,,C2,...,Cm-1 € R, no todos nulos, tales que 


ma = dl ela to ctra la) eel, 


lo que contradice la hipótesis de partida. 

<=) Supongamos que 3x0 € 1t.q.W (fi, f2,..., fm) (20) 40. 
Razonando de nuevo por reducción al absurdo, supongamos que (fi, fa,..., fm) es lineal- 
mente dependiente. Entonces, 3cC1,C2,...,Cm E R, no todos nulos, tales que 


O=c1f (2) +c2f2[2)+...Cmfm(t) Vx l, 


y derivando, 
0 = cifilo) +0 f(x) +..-Cmfml2), 


0 = aferra Da, 


para todo x E Í; esto implica que el determinante de la definición del Wronskiano es nulo 
para todo x € Í lo que contradice la hipótesis de partida. Mm 


Ejemplo 1.9 Estudiar si son linealmente independientes los siguientes conjuntos de fun- 
ciones en IR: 


(a) (pr, pa, p3) = (1, 1,1% 


=2%0. 


Our 
N 
3 


1 
W (1, P2, P3) (1) =|0 
0 


Por tanto, (1, 1,1?) es un conjunto linealmente independiente. 
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(b) Lf, f2) = Le”, e] 


e 
e 


W (fi f2) (2) = 


e 
E Apt 


dx 


=3I0 0, 


Por tanto, ([e”, e*”) es un conjunto linealmente independiente. 


(c) (br, d2, 93) = Le”, ze”, a%er) 


si ser qler 1 2 q? 
W(d,0,03) (1) = | er (x+1)e? (1 +2) xe” =e* [1141  22+2% 
e (x+2)e" (22+40+2)e* 1 242 2274242 
1 
= e*|01 2x2 |=e* (44+2-4x)= 28% 20. 
0 2 4142 


Por tanto, [e”, ze”, 12e”Y es un conjunto linealmente independiente. 


(d) (91, 92) = (sen x, cos 2) 


sen COS TL 
COST  — Sen. 


W (91, 92) (2) = 


= — sen? x — cost =-—1%0. 


Por tanto, [sen x,cos xp es un conjunto linealmente independiente. 


(e) [0,,02,03,04) = [1?—1,1,22+x,12+1) 


qe=1 
2x 
2 
0 


W (91, 0, 03, 04) (2) == 


1 22+x +1 

0 2x+1 1 =0 
0 2 0 ES 
0 0 0 


Por tanto, (1? — 1,1,1?+ x,1 +1) es un conjunto linealmente dependiente. 


Lema 1.10 Sean u;¡,us,...,un n soluciones de (Py) en el intervalo 1. Si 3 
W (u,,uz,...,Un) (19) =0, entonces W (uy, us,... 


Demostración. 


Supongamos que Jxp € Í tal que W (u;, uz, ... 


R, no todos nulos, tales que 


c1U1 (Lo) + c2ua (Lo) +- 
Derivando n — 1 veces, 
cu, (Lp) + c2u) (Lo) +: 


(n—1) 


, Un) (19) = 0, entonces existen C1,C2,... 


«++ CpUn (Lp) =0. 


+ Cpu, (tp) =0, 


cu, (20) + cura (10) +:--+ cur (To) =0: 


zo E Í tal que 
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Sea u (1) = ciu (1) + caua (1) + -*- + Cptin (1), entonces, u(1p) = u (tp) = +: = 
ulr=D (209) = 0. Además, puesto que uz, uz, ..., un son soluciones de (Pr), entonces, u es 
solución de (Py) (Proposición 1.3). Por tanto u verifica el P.V.T. asociado a (Py), con 
u(zp)) = uv (20) = +: = ub=D (2) = 0. Puesto que el P.V.I. tiene solución única por el 
Teorema 1.1, se tiene que u(1) = 0 para todo x € [, lo que implica que ([u¡,uz,..., Un) 
es un conjunto de funciones linealmente dependiente y por tanto W (u¡,uz,...,Un) (1) = 
0 VIEcl. a 


Teorema 1.11 Sean u1,u2,..., Un N soluciones de (Py) en el intervalo I. Entonces, 
(uy, uz, ...,Un | es linealmente independiente => W (u;z,uz,...,Un) (1) 40 Vx € 1. 


Demostración. 

>>) Supongamos que (uz, uz,..., un) es un conjunto de funciones linealmente indepen- 
diente. 

Razonando por el método de reducción al absurdo, supongamos que J3xp € Í tal que 
W (u¡,uz,...,Un) (10) = 0, entonces, por el Lema 1.10, W (u¡,uz,..., Un) (1) =0 Vx € 
I, lo que contradice la hipótesis inicial. 

<=) Supongamos que W (u¡,uz,...,Un)(1) 40  Vre€El. 

Utilizando la Proposición 1.8 se tiene que ([u¡, uz,..., un) es linealmente independiente. 


nm 
Teorema 1.12 Siempre existen n soluciones linealmente independientes de (Py) en 1. 


Demostración. 
Sea Tp E l y sea u;, con 1 <1< nm, la única solución del P. V.I. de la ecuación diferencial 
(Pg) con la condición inicial 


(5) Jl sij=:8t-1 
ze 2) a sijHi1 


para 0 << n, donde yo = u;. Como 


10. 0 
0 dE 
Wu ta) (e) = ra Y SiO 
Oia Y 
entonces, por el Lema 1.10 y el Teorema 1.11, (u¡,uz,..., un) es un conjunto de n soluciones 
linealmente independientes de (Py) en 1. € 
Definición 1.13 Cualquier conjunto [u,,uz,...,un y de soluciones de (Py) que sea lineal- 


mente independiente en I se dice que es un conjunto fundamental de soluciones de 
(Py) en 1. 


Teorema 1.14 Sea [u¡,uz,...,un) un conjunto fundamental de soluciones de (Py) en 1. 
Entonces, cualquier solución de (Py) se puede escribir como combinación lineal de 
U1,U2,..., Un- 
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Demostración. 
Sea u una solución de (Py), 70 E I y sean a; = ul) (2) para 0 < j < nm. Consideremos 
una función w definida como 


w(x) =c1u1 (1) + cu (1) + :>* + pun (1) Vx € l (1) 


cOn C1,C2,...,Cn € R. Buscamos los valores de estas constantes de tal forma que la función 
w verifique que 04) (27) = a. Entonces, para xy se tiene que 


w(zo) = cru (Lo) + C2ua (20) ++ + CnUn (20), 
w'(2o) = ciu (Zo) + cas (Lo) + *** + Cpu, (Lo), 
(n—1) — (n—1) (n—1) A (n—1) 
w (20) = Cu, (to) +c2ua (20) +: +epur 7 (2o), 
es decir 
uz (2p) uz (Zo) Ur, (20) C1 00 
uy (2p) u2 (Zo) Y, (To) Co | 1 
ye (20) ue (20) --* a (20) Cn An-—1 
Este sistema tiene solución única para C,,C2,...,C. ya que W (u¡,Uz,...,Un) (20) 40 
por ser (u¡,uz,...,u,) un conjunto fundamental de soluciones. 
La función w es solución de (Pr), por ser combinación lineal de (u¡, uz,..., Un), y además 


verifica que 04) (29) = a para0<j<n. 

Por otro lado, la función u también es solución de (Pr) y verifica que u% (2p) = ay para 
0<3<mn. 

Por tanto, en virtud del Teorema 1.1, u(x) = w(x) para todo El. Mu 


Observación 1.15 Del Teorema 1.14 se deduce fácilmente que no pueden existir más de 
n soluciones de (Py) linealmente independientes puesto que cualquier solución de (Py) se 
puede escribir como combinación lineal de n soluciones de (Py) linealmente independientes. 


Observación 1.16 En virtud de la Proposición 1.5 y este último resultado, se deduce que 
para determinar todas las soluciones de (P), basta con encontrar una solución up de (P) y 
un conjunto fundamental de soluciones de la ecuación diferencial homogénea asociada (Pp), 


(uz, uz,..., Un). De esta forma, cualquier solución u de (P) se puede escribir como: 
u(x) ==up (1) + cqu1 (1) + cqU (1) + --* + Cptn (1) Vx € l (2) 
con C1,C2,...,Cnh constantes. A la expresión uy (1) = C1U1 (1) + Cqua (1) +-:* + Cy Un (1) se 


le denomina solución general de la ecuación diferencial homogénea, (Py), y a la 
expresión (2), solución general de la ecuación diferencial completa, (P). 


2. Ecuaciones diferenciales lineales homogéneas con coeficientes constantes dl 


2. Ecuaciones diferenciales lineales homogéneas con coe- 
ficientes constantes 


Una vez que conocemos la estructura de la solución general de una ecuación diferen- 
cial lineal, veremos como calcular un conjunto fundamental de soluciones de una ecuación 
diferencial lineal homogénea, que como ya hemos demostrado en el Teorema 1.14 permite 
obtener la solución general de (Py). Encontrar un conjunto fundamental de soluciones de una 
ecuación diferencial lineal homogénea depende en gran medida de los coeficientes de ésta. 
Estudiaremos únicamente las ecuaciones diferenciales lineales homogéneas con coeficientes 
constantes. 

Si el orden n de la ecuación diferencial es igual a uno, tenemos la ecuación diferencial de 
primer orden 

au! (1) + bu (1) =0 


(con a,b € R constantes), o equivalentemente 


b 


ul (2) = —Zu(e) = quí) 


b 
donde y = ——. La solución general de esta ecuación diferencial es 
a 


WMul= Ke, 


con K ER. 
Para el caso n = 2, la ecuación diferencial será de la forma 


au” (1) + bu! (2) +cu(x) =0 (3) 


(con a,b,c € R constantes). Un conjunto fundamental de soluciones de esta ecuación di- 

ferencial debe tener dos soluciones linealmente independientes. Una posibilidad es ensayar 

con funciones similares a las obtenidas para n = 1. Así, podríamos pensar que la función 

u (1) =e*"” podría ser solución de (3) para algún valor A € R. 
Estudiar si esta función es solución es sencillo: 


au” (2) + bu' (1) + cu (a) =aAe*” + b4e%” + 0e%* = (a +bA +0) e%* =0. 


Puesto que e?” 4 0 para todo x € R, la única posibilidad para que u (1) sea una solución es 
que 
A A 


Por tanto, esta ecuación determina los valores de A que hacen que u (2) = e*” sea solución 


de (3). A esta ecuación se le denomina ecuación característica. La ecuación característica 
viene dada por un polinomio de grado 2 en A, 


p(A) =a A +bdA+c<c, 


que se llama polinomio característico, y que, por tanto, puede tener dos soluciones reales, 
una única solución real (de multiplicidad dos) o dos soluciones complejas conjugadas. La 
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expresión analítica de dos soluciones linealmente independientes dependerá de las soluciones 
de la ecuación característica: 


q 2 
E E E a a 


2a 
Así: 
= Sib? — 4ac > 0 => 311,42 ER, con A1 X% Az, soluciones de la ecuación característica; 
es decir, 
dla = et 
uo lx) = er 


son dos soluciones de (3). Además son linealmente independientes, pues 


A¡T Maz 


e e 
Az Maz 
Aje Ae 


Was | = (AJA o, 


Luego di es un conjunto fundamental de soluciones y por tanto la solución 


general de (Pp) es 


Aj7 A27 


u (1) =cre"” + ce 


ñ Z b : A E j 
= Sib? — 4ac=0 => 311 = e E R solución de la ecuación característica; es decir, 
a 


u (1) =e%* 


es una solución de (3). En este caso tenemos que encontrar otra solución de (3). En- 
sayamos con 


uo (1) = xe** 
y comprobamos que también verifica (3). Para ello calculamos las dos primeras derivadas 
de us: 
ala = eze" =(1+J0) e” 
un(a) = Ae" +(1+Ax) Ae" = a +2) 0? 
y por tanto, 


aus (2) + bus (2) + cuz (a) = (a(22+Wx) +b(1+Ax) + cx) e** 
= (2 (a +bA+c) + 244 + b) En 


(20 b)e* = (2 (=2) E ») e*=00 


puesto que A = —— es solución de la ecuación característica. Además, uy y us son 


a 
linealmente independientes. En efecto, 

Ar qert 
Ae” (14+A2)e%” 


(+ Ao= 47) = 420. 


1 x 
A 1+Az 


24x 


W (u,,ua) (1) = =6 
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Por tanto, Len erat es un conjunto fundamental de soluciones de (3) y la solución 
general de (Pr) es, en este caso, 


u (1) = ce" + core”. 


= Si b?— 4ac<0 => 311,42 € C conjugados que son solución de la ecuación carac- 
terística y por tanto 


> A =Q+ Bi 
Ja, BER tal que Í A 


Puesto que A y A2 son números complejos, 


w (2) = e” 
wo (1) = er? 


son funciones complejas que son solución de (3). A partir de ellas podemos buscar dos 
funciones reales de variable real que también sean solución de (3). Por la Proposición 
1.3, sabemos que cualquier combinación lineal de ww, y wa también es solución de (3). 


1 1 


u(a) = 5 (0) + 5w2 (0) = 5 (4? 4 02) 


1 
2 2 2 
_ 1 ax Bai ax —BriYy 1 ar Bai —Bxi 
= (eve + ee E (0% y 674) 
1 
= q (cos (Bx) + ¿sen (Bx) + cos (—Px) + ¿sen (—Bx)) 


= e cos (Bx) 


ua €) A 0 (x) + a (x) = > Adi _ Ets) 


AE - (ent ta A een E, en (cbr e er) 


= Ge (cos (Bx) + ¿sen (Bx) — cos (—Px) — ¿sen (—Bx)) 
=  —e* sen (Bx) = e” sen (Br) 


son también dos soluciones de (3). Además, u, y uz son dos funciones reales y son 
linealmente independientes: 
ee” cos (Bx) e”? sen (Bx) 


e” (a.cos(Bx) — PBsen (Bx)) e” (asen (Bx) + PB cos (Bx)) 
e** (asen (Bu) cos (Ba) + B cos” (Bu) — asen (Bx) cos (Bu) + PB sen” (Bu)) 
= pee + 0. 


Por tanto, [e*” cos (Bx) ,e*” sen (Bx)) es un conjunto fundamental de soluciones de (3) 
y la solución general de (Py) es 


W (uz, u2) (2) = 


u (1) = cie” cos (Bx) + cae”” sen (Bx). 
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Ejemplo 2.1 Calcular la solución general de las siguientes ecuaciones diferenciales ho- 
mogéneas de coeficientes constantes de orden 2 . 


(a) u” (1) + u' (2) — 6u (1) =0 


Calculamos primero las soluciones de la ecuación característica: 


-1+y1-4(-6) -1+5 AL=2 
2 ÑRÑ = en 1 ) 
AMHFA=6=0=>A > 3 OS 


por tanto, un conjunto fundamental de soluciones vendrá dado por las funciones: 
a A 
tala) = €, 
y de esta forma, la solución general será: 
u (1) = eye? + eye”. 
(b) y” — 4y' + 4y =0 


Empezamos calculando las raíces del polinomio característico: 


4+y16-4-4 
AM—4A+4=0>A= ; = 2 (doble), 


por tanto, las funciones 


forman un conjunto fundamental de soluciones y la solución general de esta ecuación dife- 


rencial será: 


y (1) = cre* + coxe” = (0, + aca) e”. 


(c) u” (2) +u(1) =0 
Resolvemos la ecuación característica: 


M41=0+A= kV 1 = ti > lar 
DA 


es decir, un conjunto fundamental de soluciones está formado por 


WlE) == 008%, 


ul) = senz, 
y por tanto, la solución general viene dada por: 


u (1) = cc, cos x + C2Senz. 
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(d) u” (12) — 6u' (2) + 13u (1) =0 
Calculamos las raíces del polinomio característico: 


6+y36-4-13 6+y-16 
2 o 3 


A =3+21, 


2_ == == 
*—-61+13=0=>A a 


=3+00>( 


así, un conjunto fundamental de soluciones viene dado por 
u (2) = e*cos(2x), 
uo (1) = e*sen(21), 
y por tanto la solución general viene dada por: 
u (1) = cre* cos (27) + cae” sen (2x) . 


En el caso general de ecuaciones de orden n el procedimiento que se sigue es el mis- 
mo. En este caso tendríamos que dada una ecuación diferencial homogénea con coeficientes 
constantes, 

au” (0) +---+agu" (2) + aju' (2) + ayu (a) =0, 


se probaría, de forma análoga al caso de orden dos, que el polinomio característico es 


P(A) =0,M + +02 + a1A + 00, 


y a partir de sus raíces podremos construir un conjunto fundamental de soluciones: 


= Por cada raíz real simple A del polinomio característico, la función v (1) = e*” forma 
parte del conjunto fundamental de soluciones. 


= Por cada raíz real A de multiplicidad k del polinomio característico, las funciones 
v¡(u) =x1e%*, con 0< ¡< k, forman parte del conjunto fundamental de soluciones. 


= Por cada raíz compleja y su conjugado, A = a. = 61, del polinomio característico, 
las funciones v1 (1) = e cos (Bx) y va (1) = e” sen (Bx) forman parte del conjunto 
fundamental de soluciones. 


= Por cada par de raíces complejas y su conjugado, A = +63, de multiplicidad k del poli- 
nomio característico, las funciones 0,1 (1) = a1e"* cos (Bx) y vja (1) = 21e** sen (Bu), 
con 0 < j < k, forman parte del conjunto fundamental de soluciones. 


Ejemplo 2.2 Calcular la solución general de la ecuación diferencial homogénea con coefi- 


cientes constantes 
u” (2) + 3u” (2) — 4u(x) =0. 


En este caso la ecuación característica es 243124 = 0 y tiene por solución Ay = 1 (sim- 
ple) y Az = —2 (doble). Por tanto, un conjunto fundamental de soluciones es [e*,e7%, ey, 
y la solución general es: 


u (2) = cre? + ce Y + cgre ”. 
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Ejercicios propuestos 
Calcular la solución general de las siguientes ecuaciones diferenciales homogéneas de 


coeficientes constantes: 


1. 3y" — 19y” + 36y' — 10y =0 Solución: y (1) = c,e3% + cseó cos x + eze” sen z. 


2. u” (2) +24" (2) +u(x) =0 Solución: u (1) = €, COS Y + C2SEN LT + (31 COS LT 
HF CATGEN TZ. 
3. u' (2) —10u (1) +29u(1) =0 Solución: u(x=) = ce "cos(2x) + cae*” sen(27). 


4. u"(x) — 3u (1) + 2u(1) =0 Solución: u(x) = ce % + cae” + czxe”. 


3. Método de los coeficientes indeterminados 


Recordemos que para resolver una ecuación diferencial lineal de orden mayor que uno 
necesitamos encontrar un conjunto fundamental de soluciones de la ecuación diferencial ho- 
mogénea asociada y una solución particular de la no homogénea. En la Sección 2 hemos visto 
cómo calcular un conjunto fundamental de soluciones de (Py). A continuación describiremos 
un método que nos servirá para calcular una solución particular de la completa. 

En general, encontrar una solución particular de una ecuación diferencial es un trabajo 
muy complicado; sin embargo, si nos fijamos en la estructura de la ecuación diferencial que 
estamos estudiando, en muchos casos será posible ensayar con una función similar al término 
independiente, b(x1), que nos lleve a una solución particular de la ecuación diferencial no 
homogénea. 

Veamos algunos ejemplos. 


Ejemplo 3.1 Calcular la solución general de la siguiente ecuación diferencial: 
ul” (1) + 4ul (2) — 2u (1) = 21? — 3x +6. 
Sabemos que la solución general de esta ecuación diferencial será: 
u (1) =up (1) + c1uy (2) + cqua (1), 


donde f(u,,uz) es un conjunto fundamental de soluciones de la ecuación diferencial ho- 
mogénea asociada, (Py). Calculamos u, y us; para ello, resolvemos la ecuación 


uy (1) + 4uy, (1) — 2ug (1) =0. 
Por tanto, p(A) = 2 +4A — 2, y sus raíces vendrán dadas por 


Mi 
z ; a 


HA) =0=>A4 
es decir, 
uy (1) = ce 2+ 6) + eje (2 V6)2. 


Únicamente faltaría calcular up (1). Para ello, puesto que el lado derecho de la ecuación 
es un polinomio de grado 2, resulta lógico pensar que una solución up pueda ser un polinomio 
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del mismo grado. Por tanto, ensayamos con un polinomio de grado 2 con unos coeficientes a 
determinar, es decir, 
up (1) = Ax* + Bx+C. 


Comprobamos si existe algún valor para los coeficientes A, B,C' de tal manera que up sea 
solución. Para ello necesitamos calcular las dos primeras derivadas de esta función y sustituir 
en la ecuación diferencial: 


uplx) = 241+B, 
ue (e) == DA; 


y por tanto, 


up (1) + 4u, (1) — 2up (2) 


(24) + 4(24x + B) -2(4x? + Bx+C) 
= 241? + (84 -2B)x + (24 +4B -2C). 


Este polinomio así calculado y b(x) son iguales, 
—2Ax? + (8A-2B)1 + (24+4B-2C) =22*— 32 +6, 


de modo que, identificando coeficientes, llegamos al siguiente sistema lineal de ecuaciones: 


—2A ==) 2 
ga -—2B = —3 
24 +4B -2C = 6 


cuya solución es A=-—1, B =-—5/2 y C = —9. Por tanto, 
up (1) =-xw*-—Zxr-—9 
2 
y la solución general de la ecuación diferencial completa será 
5 
u (1) = —x? — ¿a 9+ eje (-2+v6)8 + eye 2-V8)=, 


Si nos fijamos en la manera de calcular la función up en este ejemplo, la idea ha sido 
sencilla: buscamos una función que al sustituir en la ecuación diferencial (derivando, suman- 
do, etc...) nos dé un polinomio. Obviamente, una posibilidad (aunque no es la única) es un 
polinomio del mismo grado. Como veremos a continuación, si el lado derecho de la ecuación 
b(x) es una función trigonométrica, tendremos que ensayar con una función trigonométrica; 
si fuera una exponencial, ensayaríamos con una función exponencial, etc... La función con la 
que ensayemos siempre dependerá de unos coeficientes que tenemos que determinar. Á este 
método se le conoce como método de los coeficientes indeterminados. 


Observación 3.2 Con este método hemos calculado una solución particular de la ecuación 
diferencial. Como veremos en la siguiente sección, existe otro método para el cálculo de solu- 
ciones particulares que podría llevarnos a distintas soluciones. Sim embargo, pensemos que 
estamos expresando la solución general de la ecuación diferencial de una forma concreta que 
depende de una solución particular y un conjunto fundamental de soluciones. Por supuesto, 
esta representación no es única. 
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Ejemplo 3.3 Calcular la solución general de la ecuación diferencial 
u” (2) —u' (2) + u (2) = 2sen (32). 
Igual que antes, empezamos resolviendo la ecuación diferencial homogénea, asociada 
uz (1) — Uy (2) + un (2) =0, 


cuya ecuación característica es 


MA=A+1=0=>A= 


1+w1-4 Ye. 
= 2. 
) q O 


Por tanto 


e 

[50] 

O 

Il 

y 

NN 

[02] 

o) 

S 
AS 
al vl5 

2 
xs 


es decir, 


A v3 E v3 
uy (1) = c,e? cos 3% + ce? sen ls 


Buscamos ahora una solución particular. En este caso, puesto que b(x) = 2sen (3x1), resulta 
lógico pensar que la solución particular puede ser una combinación lineal de las funciones 
sen (31) y cos (31). Así, ensayamos con: 


up(x) = Asen(3x) + Bcos (31) 
up(x) = 3Acos (3x1) — 3B sen (3x1) 
up(x) = —9Asen(3x) — 9IBcos (31), 


y sustituyendo en la ecuación diferencial, se tiene que 
up (2) — up (2) + up (1) = (94 +3B + A) sen (32) + (-9B — 34 + B) cos (32). 


Por tanto, identificando los coeficientes de esta última expresión con los de b(x) = 
2sen(3x), nos queda el siguiente sistema lineal de ecuaciones en A y B: 


=8A +3B = 2 
=3A —8B = 0 


cuya solución es A=-—= y B = S La solución general de la ecuación diferencial completa 
es 


16 6 2 3 2 3 
uu) = 73 Sen (31) + 73 008 (31) + cie? cos ES + ce? sen (E) : 
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Ejemplo 3.4 Calcular la solución general de la ecuación diferencial 
y” —2y' — 3y = 4x0 —5+6zxe”. 
Empezamos resolviendo la ecuación diferencial homogénea asociada, cuya ecuación car- 


acterística es 
2+ v4+12 = 
Momoa0> 1 EE a 
2 Ti, 


y por tanto 


y (a) = cre** + eze”. 


Puesto que b(1) es la suma de un polinomio de grado 1 y del término ze”, buscamos una 
solución particular ensayando en este caso con 


yr (1) = Ar + B+ Ce” + Dxe” 
De aquí, 


yr (a) = A+2Ce* + De* + 2Dxe* = A+ (20 + D)e* + 2Dxe*, 
yo(x) = 2(2C+.D)e* +2De* + 4Dzxe* = (40 + 4D) e* + 4Drxe*, 


y sustituyendo en la ecuación diferencial, 
yr (2) —2yp (1) - 3yp (4) = =34x +(=3B -24) + (-30 +2D)e* — 3Dxe*”, 


que, como en los ejemplos anteriores, nos permite obtener para los coeficientes 4, B,C, D el 
sistema de ecuaciones 


—3A = 4 
24 -—3B = -—5 
3C +2D = 0 
—3D 6 
cuya solución es A = =3, B = 2, c= =5 y D = —2. De esta forma, la solución general de 
la ecuación diferencial es 
4 23 dl 
ulal= 37 + EN dl A E A 


Ejemplo 3.5 Calcular la solución general de la ecuación diferencial 
"” ! T 
y —0y +4y = 8e”*. 

La ecuación característica de la ecuación diferencial homogénea es 


| A=4 


2_ NE => 
A=5BA+F4=0=>A > 3 e = 


Por tanto, la solución general de la ecuación diferencial homogénea es 


Yu (1) = c,e*” + eze”. 
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Calculamos ahora una solución particular. Para ello, si observamos el lado derecho de la 
ecuación, la primera idea es ensayar con 


yp (2) = Ae” = yp (1) = yp (1) 
que si sustituimos en la ecuación diferencial nos queda: 
yr (1) — 5yp (1) + 4yp (12) = (A-54+44)e” =0. 


En este caso no existe para ningún valor de A una solución de la ecuación diferencial no 
homogénea. En realidad, podíamos habernos dado cuenta de que yp (1) = Ae” no es válida 
para ensayar, puesto que es una solución de la homogénea. 

Buscamos otra posible solución de la ecuación diferencial no homogénea. Seguimos en- 
sayando con funciones que dependan de e”, por ejemplo, 


yp (1) = Axe” 
para la que se tiene que 


yp(x) = Ae” + Axe”, 
y, (1) = Ae” +A4e” + Axe” = 24e” + Aze”. 


Al sustituir en la ecuación diferencial, 
yr (1) — 5yp (1) + 4yp (1) = (24—54)e” +(4-—54+44) xe” =-34e”, 


e igualar a b(x) =8e”, obtenemos A =—3. Por tanto, la solución general es 


8 
us sd bae? Hat; 


4. Método de Lagrange o de variación de las constantes 


Existen ecuaciones diferenciales lineales no homogéneas para las que ensayar con alguna 
función similar a la que aparece en el término independiente, b(x), no proporciona ningún 
resultado. Por ejemplo, para la ecuación diferencial 


u” (1) + u (1) = cosec z, 


ensayar con combinaciones lineales de funciones trigonométricas de tipo sen x, cos x o cosec x 
no permite encontrar up(x). 

En estos casos tendremos que recurrir a otro método que nos permita determinar algu- 
na solución particular de la ecuación diferencial completa: el método de Lagrange, o de 
variación de las constantes, que es una generalización del método utilizado para resolver 
ecuaciones diferenciales lineales de primer orden. 

En general, dada la ecuación lineal de segundo orden 


as (1) u” (2) + ay (xju! (2) + ay (12) u(2) =b (2) 


4. Método de Lagrange o de variación de las constantes 1% 


definida en el intervalo / € IR, su ecuación homogénea asociada tiene dos soluciones li- 
nealmente independientes que denotaremos por u, y us. La solución general de la ecuación 
diferencial homogénea es por tanto 


uy (1) = c1u1 (1) + caua (1). 


Igual que en el caso de las ecuaciones difereciales lineales de primer orden, para calcular 
la solución general de la ecuación no homogénea, sustituimos las constantes C¡ y C2 por 
las funciones c1 (1) y cz (1) respectivamente. Nuestro objetivo será ahora calcular dichas 
funciones c1 (1) y cz (1) para las que 


up (1) =C1 (12) u1 (1) + ca (1) ua (x) 
sea solución de la ecuación diferencial no homogénea. Por tanto, derivando 
up (a) = 0 (1) u (2) + cz (2) ua (2) + cr (2) u (1) + cz (1) uz (2) 
up (1) (6 (2) 1 (2) + 6 (1) us (1) +01 (2) us (1) + e) (2) us (1) 
+c1 (2) uy (2) + c2 (2) us (1) 


y sustituyendo en la ecuación diferencial 


b(x) = az(x)u” (2) +a,(vJu' (2) + ay (2) u (1) 
= a (2) [(c, (2) 41 (2) +0) (0) ue (1)) +0 (2) us (2) + 6, (2) us (a) 
+42 (2) [c1 (20) u (2) + cz (2) ue (2)] 
+41 (1) [c, (2) us (2) + ca (2) ua (2) + 01 (2) ul (2) + ca (2) us (2) 
+0 (1) [c1 (1) ur (1) + ca (2) ue (2)] 
= C1 (2) la, (2) y (2) + ar (2) u, (2) + an (2) us (2) 


donde hemos utilizado que u¡ y us son soluciones de la ecuación diferencial homogénea 
asociada. Si calculamos todas las funciones cy (1) y ca (1) que verifican 


b(2) = ar (o) [(c, (2) uu (2) + ca (2) ue (2) + (2) u, (2) + 0, (1) uz (2) 
+41 (2) [c, (2) u (1) + c, (2) uz (1)] 


tendremos la solución general de la ecuación diferencial completa. Sin embargo, sabemos que 
para calcular la solución general basta con obtener una solución particular de la ecuación 
diferencial no homogénea y un conjunto fundamental de soluciones de la homogénea (que ya 
conocemos, (u,, uz)). De esta manera, es suficiente calcular una función c,(1) y otra ca(x), 
que verifiquen esta última ecuación. Por tanto, si c(1) y ca(1) verifican 


cy (2) us (2) +0) (u)uo (1) = 0 


b(x) (4) 


cy (0) uy (2) +02 (2) us (1) 
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tendremos que up (1) =C1 (1) u1 (1) + ca (1) ua (1) es una solución particular de la ecuación 
diferencial completa. 

Observemos que para calcular c¡ (1) y ca(1), puesto que uy y uz son funciones conocidas, 
tenemos, para cada x E I, un sistema lineal de dos ecuaciones con dos incógnitas bien 
definido, pues as (1) 4 0 para todo x € [. Además, dado que 


W (ur, uz) (1) = dd (1) us 


E | + para todo x € I, 


al ser (uz, uz) un conjunto linealmente independiente, se tiene que el sistema de ecuaciones 
(4) es compatible determinado. Resolviendo por el método de Cramer: 


de us (1) u (1) bla 
— ue (e) uy (a) E 
NO BIOS, + ¿us aa () |. 65) 


con lo que 


Observación 4.1 En el caso de ecuaciones diferenciales de orden superior el procedimiento 
es análogo. Así, para la ecuación diferencial 


az (2) u” (2) + a2 (2) u” (2) + ar(xju' (2) +a0(2)u(x)=b(x), 
una solución particular vendrá dada por 
up (1) =C, (2) u1 (1) + c2 (12) ua (1) +03 (1) uz (1), 


siendo [u,,uz,uzy un conjunto fundamental de soluciones de la ecuación diferencial ho- 
mogénea y Cc (1),c2(1) y c3 (1) funciones que verifican 


4 (0) () + (0)ue(a) + dleju(e) = 0 
alu le + dll) + ál)ul) = 0 
40 () + (0) + alu = Ll 


Ejemplo 4.2 Calcular la solución general de la ecuación diferencial 
u” (2) + u (1) = cosec z. 


Como en los ejercicios anteriores, empezamos resolviendo la ecuación diferencial ho- 
mogénea asociada, uy (1) + uy (1) = 0, cuya ecuación característica es: 


X2+1=0>A=wi, 
y por tanto 


u (1) = cosz, 


us(1) = senz. 
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Entonces, la solución de la ecuación diferencial homogénea asociada es 
Uy (1) =C, COST + C2SEN 2. 
Utilizando el método de Lagrange, buscamos una solución particular up tal que 
up (1) =Cc, (1) cos x + ca (1) sen z. 


Las funciones c,(1x) y ca(x) tienen que verificar el sistema (4), que en este caso, es: 


E (E) cos E + c)(1)senz 


0 
(1) (—senxz) + cs(x)cost = cosecz 1 


Aplicando el método de Cramer obtenemos la solución que viene dada por (5). Por tanto, 
calculamos primero el Wronskiano de Lu;, us]: 


COS XT sen Y 


= cos? x + sen? x = 1, 
— sen COST 


W (u1,u2) (1) = | 


luego, 
0 sen x 
'(a) | COSECT COST : 
ea = — sen 1 cosec x= — 
1 
W (u;¡, us) (2) 
ya 
COS Y 0 
; —senx  cosecz 
es (== = COS T COSeC Y = COt$8 1. 


W (uz, u2) (2) 


Obtenemos sus primitivas, 


e (1) = f san =-x 


ola= foot xdx = log (sen 1). 
De esta manera, la solución general de la ecuación diferencial completa es 


u (1) =—ucos x + log (sen 1) sen x + C1 Cos 1 + C2¿SEnz. 


Ejemplo 4.3 Calcular la solución general de la ecuación diferencial 
w” (2) +u (1) =tgz. 


Empezamos calculando la solución de problema homogéneo asociado: 


U= AGA ZA M1) >< A, 


> 
[o] 

! 

| 
5 
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por tanto, un conjunto fundamental de soluciones está formado por las funciones u¡ (1) = 1, 
uz(x) =cosx y uz(x) = sen x, es decir, 


Uy (1) =C1 + C2C08 7 + C3Ssenz. 
Utilizando el método de Lagrange, buscamos una solución particular de la forma 
up (1) = c1 (1) + cz (1) cos x + c3 (1) senz 


y calculamos c, (1), ca (1) y cz (1) que verifican el sistema 


(1) + c(x)cosx + (4 (x)senzx =00 
co (x)(=senx) + ((1)cosz 0) 
Cc (1) (=cos1) + «G(r)(—senz) = tgz 


que resolvemos mediante el método de Cramer. Para ello, obtenemos primero 


l  cosz sen z 
W (u,,uz,uz) (1) =| 0 —senx cos |=1 
O —cosr —senz 


y a continuación 


0 COS Y sen x 
O —senx  cosz 
da ter —COost —Senz 

W (ur, uz, uz) (2) 


=tg8 1, 


1.0 sen z 
0 0 COS YT 
O tex —senz 


ela = —COSIÍE X= — SON L, 
2) = 72 (a) á 


l  cosz 0 
O —senx 0 
O —cosx tez 


W (u;, ua, uz) (1) 


Calculamos ahora las primitivas de estas funciones: 


frena = —log (cos 2), 


= —senztgx. 


c1 (1) 


ca (1) f- sen xdx = cos 1, 


2 2 
sen? x cosx—1 
| asenotgndr= | - dx ¡Zu 
COS Y COS 


1 1 
J(es=- ) do =wenw= | dx 
COS YT COS YT 


c3 (1) 
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Para calcular esta última primitiva, multiplicamos y dividimos por cos x y hacemos el cambio 
de variable t = sen x, con lo que tenemos: 


] > d ] . d ] z dt ] 3 + > dt 
a cos xd = | ——dt = 
COS XI cos? x 1-2? 1+t 1-t¿ 


1 | il Naas 
tac 


log 


2 a e O: 1-22 
1 1 4 1 

= las psa) A 
2 cos? x COS y 


de modo que, 
c3 (1) = sen x — log (sec x + tg 1). 


Por tanto, la solución general de la ecuación diferencial es 


u (2) =—log (cos 1) + cos? x + (senx — log(sec x + tg 2)) sen x + c1 + ca cos x + cz sen 2. 


Ejercicios propuestos 


1. u"(x) — 6u'(x) + 9u(x) = 62? + 2 
2 8 2 
Solución: u(1) =512%+-x +2 +c/e% + coue?”. 


3 9 3 
2. y" (2) 3u (1) a (2) = 14 sen(21) — 18 cos(21) 


Solución: u (1) = 2sen(2x) + 3cos(2x) + c,e* + cge””. 


3. u"(1) +4 (1) = xcosx 


1 1 1 
Solución: u(1) = a sen T +C08s 2 + zo sen x — ge cos T+C1+ cae" 


4. u"(x) + u(x) = 41 +2 


1 3 1 3 
Solución: u(1) = 4x4+2+c,e *+c2e2” cos (2) +c3e2* sen (E) 
5. u' (2) — 3u (2) + 2u(x) = e” sen(e”) 


Solución: u(x) = (e? cos(e”) — sen(e*)) e? — cos(e*)e?” + ce” + ege?”. 


5. Ejercicios resueltos con Maple 
Ejercicio 1. Resolver el siguiente problema de valor inicial: 


u" (2) + 4u (2) = cos? x 
u(0) =0 
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> restart:with(plots): 
La ecuación que tenemos que resolver es: 


> ecuacion:=diff(y(x),x,x,x)+4*diff (y(x),x)=(cos(x))72; 
d? d 
Py (a) +4 y (2) = (cos(1))? 


Buscamos, en primer lugar, la solución general de la ecuación homogénea asociada; para ello 
resolvemos la ecuación característica: 


> EC:=173+4x*1=0; 


ecuacion := 


EC :=P'4+41=0 
> solve(1EC+,11)); 
=D = 27 EL == 20 
A la vista de los valores de l, un conjunto fundamental de soluciones de la ecuación diferencial 
homogénea está formado por las funciones u1(x), u2(1) y u3(x): 


> ul:=x->1; 
ul = mel 

> U42:=x->cos(2x*x); 

u2 := 1H cos(21) 
> u3:=x->sin(2x*x); 

us := hosin(21) 

y, por tanto, su solución general se escribe como 
> ugH:=x->C1xu1 (x)+C2*u2 (x)+C3xu3(x) ; 
ugH :=x> Clul(x) + C2u2(1) + C3u3(z) 
> ugH(x); 
C1 + C2 cos(21) + C3 sin (21) 


Para encontrar una solución particular de la ecuación completa aplicaremos el método 
de variación de las constantes. Supondremos, entonces, que existen funciones c1(1), c2(x) 
y c3(x) tales que upC(x) = cl(x)ul(x) + c2(1)u2(x) + c3(x)u3(x) sea la solución parti- 
cular que vamos buscando. Estas funciones se obtendrán resolviendo el siguiente sistema 
de ecuaciones, donde clp, c2p y c3p representan las derivadas de cl(x), c2(x) y c3(x), 
respectivamente. 


> ecul:=c1px*ul(x)+c2pxu2(x)+c3pru3(x)=0; 
ecul := clp + c2p cos(2 1) + c3p sin(2x) =0 
> ecu2:=c1px*diff (u1(x),x)+c2px*xdiff (u2(x) ,x)+c3p*diff (u3(x),x)=0; 
ecu2 := —2c2p sin(2x) +2 c3p cos(21x)=0 
> ecu3:=c1px*diff (u1(x),x,x)+c2p*xdiff (u2(x),x,x)+ 
c3p*diff (u3(x),x,x)=(cos(x))72; ecug := —4c2p cos(2x) — 4c3p sin(2x): 
> solve(fecul,ecu2,ecu3)+,tfc1p,c2p,c3p») ; 
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(cos (1))?  cos(2x) 
(cos (2 1))* + (sin (2 1)) 


3, C3p = —1/4 


(es = 1/4 (cos (a))*, c2p = —1/4 oi nO >) 
> simplify(%); 


[cp = 1/4 (cos (1), c2p = —1/4 (cos (1))* (2 (cos (a) — 1), c3p = —1/2 sin (x) (cos (137) 


Ahora integramos clp, c2p y c3p para obtener las funciones cl(x), c2(1) y c3(x). 
> Int((1/4)*cos(x)72,x)=int((1/4)*rco8(x)72,x); 
fr (cos (1))? de = 1/8 sin (%) cos (1) +1/8x% 
> c1:=x->(1/8)*rsin(x)*cos(x)+1/8xx; 
cl := xH 1/8 sin (1) cos (1) + 1/8x 
> Int(-(1/4)*rco08(x)72*(2xco0s(x)72-1),x)=int(-(1/4)*cos(x)72x*(2xco05s(x)72-1),x); 
fan (cos (a))* (2 (cos (a) — 1) de = —1/8 sin (2) (cos (1))? — 1/16 sin (2) cos (a) — 1/16x 
>  c2:=x->-(1/8)*sin(x)*cos(x)73-(1/16)*sin(x)x*cos(x)-(1/16)xx; 
c2 := 25 —1/8 sin (1) (cos (1))? — 1/16 sin (1) cos (a) — 1/16x 
> Int(-(1/2)*rsin(x)*co8(x)73,x)=int(-(1/2)*sin(x)x*rcos(x)73,x); 
fan sin (1) (cos (1))* de = 1/8 (cos (a))* 
> c3:=x->(1/8)x*cos(x)74; 
c2 := a > 1/8 (cos (%)) 
La solución particular de la ecuación diferencial completa será, entonces, 
> upC:=x->c1(x)x*ul (x)+c2(x) *u2(x)+c3(x)x*u3(x) ; 
upC := x% >cl(x)jul(x) + c2(1)u2(x) + c3(x)u3(x) 
> pc): 
1/8 sin (1) cos (1) + 1/81 + (-1/8 sin (2) (cos (1))? — 1/16 sin (1) cos (1) — 1/16 2) cos (21) + 
1/8 (cos (1))* sin (21) 
y la solución general de la ecuación diferencial completa vendrá dada por 
> ugC:=x->upC(x)+ugH(x) ; 
ugO := x >upC(x) +ugH (1) 
>. MBE) 
1/8 sin (2) cos (2) + 1/81 + (-1/8 sin (2) (cos (1))% — 1/16 sin (1) cos (a) — 1/16 1) cos (21) + 
1/8 (cos (1))*sin(2x) + C1 + C2 cos(2x) + C3 sin (21) 
Podemos comprobar que, en efecto, esta expresión obtenida para ugC(x) verifica la 
ecuación diferencial propuesta: 


> simplify(diff(ugC(x),x,x,x)+4*diff (ugC(x),x)-(cos(x))72); 
0 
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Imponiendo las condiciones iniciales, obtenemos los parámetros C1, C2 y C3 que deter- 
minan la solución particular del problema de valor inicial. 


> simplify(diff (ugC(x),x)); 


1/4 sin (1) cos (2) x — 4 C2 sin (1) cos (1) — 203 + 1/4 (cos (1)? +4 (cos (1))? 03 
>  dugC:=x->(1/4)*sin(x)*cos (x)*x-4*C2*sin(x)*cos (x) -2*C3+ 
(1/4) *cos(x)72+4x*co08(x)72x*C3; 


dugC := x > 1/4 sin (1) cos (1) a + sin (1) cos (a) — 203 + 1/4 (cos (2))7 + 4 (cos (1))? 03 
> simplify(diff(ugC(x),x,x)); 


—8 CS sin (1) cos (2) — 1/4 sin (1) cos (1) — 1/4x +4 C2 + 1/2 (cos (1))?x—8 (cos (1))? C2 


>  d2ugC:=x->-8*C3x*sin(x)x*xcos (x) -(1/4)*sin(x)x*xc0s (x) - (1/4) +*x+4*C2+ 
(1/2)*c0s (x)72*x-8*c08 (x)72*C2; 


d2ugC := 2H -8.C3 sin (1) cos (1) — 1/4 sin (1) cos (a) — 1/42 — 1+1/2 (cos (2))? x +2 (cos (a) )* 


> solve(tugC(0)=0,dugC(0)=1,d2ugC(0)=1F,1C1,C2,C37F); 


[C1 =1/4, 02 =-1/4, C3 =3/8) 


Sustituimos en la expresión de ugC(x) el valor de las constantes C1,C2 y C3 obtenidas 
> subs(1C1=1/4,C2 =-1/4,C3=3/8),ugC(x)); 


1/8 sin (2) cos (1) + 1/84 + (-1/8 sin (2) (cos (1))% — 1/16 sin (1) cos (a) — 1/16 2) cos (21) + 
1/8 (cos (1) sin (24) + 1/4 —1/4 cos (24) + 3/8 sin (2 1%) 
> simplify(%); 


6 sin (1) cos (a) — 1/8 (cos (1))x — 1/2 (cos (a:))? + 3/16x + 1/2 
y llegamos a la solución del problema de valor inicial: la función u_pui(x): 


> u_pvi:=x->(15/16)x*sin(x)x*cos(x)-(1/8)x*cos(x)72*x- 
(1/2)xc08(x)72+(3/16)*x+1/2; 


U_ pul i= 21 - sin (2) cos (1) — 1/8 (cos (a))%x — 1/2 (cos (a))? +3/16x + 1/2 
cuya representacion gráfica es 
> sol_pvi:=plot(u_pvi(x),x=-5..5,y=-2..2,color=black,thickness=3): 
> display(sol_pvi); 
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Para finalizar el ejercicio, representamos gráficamente la solución particular obtenida 
junto con algunas curvas de la familia ugC(x): 


> ugC:=x->upC(x)+ugH(x): 
> C1:=1/4:C2:=-1/4:familia:=seq(plot(ugC(x), x= -5..5,y=-6... 6, 
color=black,thickness=1), C3=-4..4): 


> display(familia,sol_pvi); 
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Ejercicio 2. Resolver el siguiente problema de valor inicial: 


u"(2) +74 (2) + 16u(2) + 12u(7) = 124? — 22 + 3e** 
u(0) = 1 

u'(0) = -3 

u”(0) =0 

> restart:with(plots): 


Definimos la ecuación diferencial que tenemos que resolver: 
> 


ecuacion: =diff (u(x),x,x,x)+7x*diff (u(x),x,x)+16x*diff (u(x),x)+12x*u(x)= 
12xx73-2*x+3xexp(-3x*x) ; 
d? 
ecuacion 


u(x) +12u (1) =123 -217+3e ** 
dx 
y resolvemos la ecuación característica para obtener la solución de la ecuación diferencial 
homogénea asociada: 
> solve(1173+7*x172+16x1+12=07),11$ ; 


> 


[1 =-—3), [i =-2), [i = -2) 


de modo que, un conjunto fundamental de soluciones de la ecuación diferencial homogénea 
estará formado por las funciones ul(1),u2(1) y u3(x) dadas por 
ul:=x->exp(-3x*x); 


u2: 


ao=eec% 
=x->exp(-2xx) ; 
4 i=axHe* 
> U3:=x->x*rexp(-2*x); 
ud := xn xe ?** 
y su solución general vendrá dada por 
> 


ugH:=x->C1xu1 (x)+C2xu2 (x) +C3*u3(x) ; 
SS 


ugH := 25 Clul(x) + C2u2(1) + C3 u3(x) 
ugH(x) ; 


A E 


Para obtener una solución particular de la ecuación completa aplicamos el método de 
> 


los coeficientes indeterminados; puesto que b(x) = 121? — 2x + 3e% y al ser e una 
solución particular de la homogénea, podemos ensayar con 


u_particular:=x->Ax*x73+Bxx"2+Cxx+D+Ex*xx*exp (-3xx) ; 
u_ particular := xo Ax? 


coeficientes A, B, C, D, y E. 


BORDE 
Sustituimos esta función en la ecuación diferencial completa para poder determinar los 
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> eval(ecuacion,u(x)=u_particular(x)); 
6A+ Ee*? 4 42 Ax + 14B + 48 Ax? + 32 Bx + 16C +12 Ax? + 12Bx*+12Cx+12D= 
10 DAS 


e, identificando coeficientes, llegamos al siguiente sistema de ecuaciones en 4, B, C, DyE: 


> solve(112x*A=12,48xA+12x*B=0,42x*A+32*B+12*C=-2,6x*xA+14*B+16x*xC+12*D=0,E=37, 


1A,B,C,D,E)); 
31 
(A 1,B + 4,0 = TD 1 O =3) 


Por tanto, una solución particular de la completa viene dada por 
> subs(1A = 1, B= -4, C =7,D = -31/6, E = 3),u_particular(x)); 


a—4+ 77 _ +3xe ** 
> UPpC:=x->x7"3-4*x"2+7*x-31/6+3*xx*exp(-3*x); 
up :i= xo 1040 + 70 _ Loge + 
y la solución general de la completa será, entonces, 
> UugC:=x->ugH(x) +upC(x); 
ugO := x>uygH (1) +90 (1) 
> uc: 
CLE*t* ECTS PBOIE Rd ao _ EEES E IN 
Podemos comprobar que, en efecto, esta expresión obtenida para ugC(x) verifica la 
ecuación diferencial propuesta: 


> diff(ugC(x),x,x,x)+7*diff (ugC(x),x,x)+16x*diff (ugC(x),x)+12*ugC(x) 


-12x*x73+2*x-3x*exp(-3*x) ; 
0 


Para resolver el P. V.I. imponemos las condiciones iniciales dadas en el enunciado, 


> diff(ugC(x),x); 
301" =26%8 0? =2031 **Lt30 87 ABS A? 


> dugC:=x->-3*C1*exp(-3*x) -2*C2xexp(-2*x)+C3x*exp (-2*x) -2*C3*xx*exp (-2*x)+ 


3*x72-8*x+7+3x*exp(-3*x) -9*xrexp(-3*x) ; 
due ¿=015=3 CL **=2.026 "LL 038072040 bd 8d TEE A e 


> diff(ugC(x),x,x); 
A E OS A AS E E E A A 


>  d2ugC:=x->9*C1*exp(-3*x)+4*C2x*exp(-2*x) -4*C3*exp (-2*x) +4*C3xxx*exp (-2*x) 


+6*x-8-18xexp(-3*x)+27*x*rexp(-3*x) ; 
d24gO ES e UCLE PE PACIEN 
> solve(tugC(0)=1,dugC(0)=-3,d2ugC(0)=07+,4C1,C2,C3H; 
(C1 = -4/3, C2 = 15/2, C3 = -2) 
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Sustituimos el valor de las constantes C1,C2 y C3 obtenidas en la expresión de ugC(x) 
> sgubs(1C1 = -4/3, C2 = 15/2, C3 = -2),ugC(x)); 


31 
—4/30 ** + 15/20 **-—21e +2 434717 ra +3 2x0 ** 
y llegamos a la solución del problema de valor inicial: la función u_ pui(x): 


> u_pvi:=x->-(4/3)*rexp(-3*x)+(15/2)*exp(-2*x) -2*x*exp (-2*x) 
+x73-4*x72+7*x-31/6+3*x*exp (-3*x) ; 


31 
4 pori= > -d/30 415/20? =2x0 “tad Tr 6 +3xe ** 
Su representación gráfica es 


> sol_pvi:=plot(u_pvi(x),x=-1..3,y=-4..2,color=black,thickness=3): 
> display(sol_pvi); 


Finalizamos el ejercicio representando gráficamente esta solución junto con algunas de 
las curvas solución de la ecuación completa. 


> ugC:=x->ugH(x)+upC(x): 
> C1:=-4/3:C2:=15/2:familia:=seqíplot(ugC(x), x= -2.. 4, y = -6 .. 6, 
color=black,numpoints=500), C3=-5..5): 


> display(familia,sol_pvi); 


5. Ejercicios resueltos con Maple 


Ejercicios propuestos 


1. Resolver el siguiente problema de valor inicial: 


uv" (2) +4" (2) + 4u (2) + 4u(x) = 81? — 24? + 1-— 1l0senx 


u(0) =—9 
uw'(0)=0 
u”(0) =-—11 
il 29 6 9) 9) 1 
Solución: u(x) = =50*=2p cos 21 — 75 Sen 2x— g $en y g osa 22% — 


2. Resolver el siguiente problema de valor inicial: 


u"(x) + u (2) = te(x) sec(x) 


ato =1 
au (0)= 
a (=1 


Solución: u(1) = usen x + cos xlog(cos 1) + 1 
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